Polynomials involving the floor function 

Inger Johanne Haland 1 and Donald E. Knuth 2 

Abstract. Some identities are presented that generalize the formula 
x 3 = 3x [x [x\\ - 3 [x\ [x [x\ J + L^J 3 + 3{x} {x [x\ } + {x} 3 
to a representation of the product xqx± . . . x n -\. 

1. Introduction. Let \_x\ be the greatest integer less than or equal to x, and let {x} = 
x — L^J be the fractional part of x. The purpose of this note is to show how the formulas 

xy = [x\y + x[y\ - [x\ [y\ + {x} {y} (1.1) 

and 

xyz = x[y[z\\ +y[z[x\\ + z[x[y\\ 

~ l x \ [yl z \\ ~ lv\ [ Z V X \\ - z[x[y\\ 

+ L^J lv\ L^J 

+ {%} {y [z\ } + {y} {z [x\ } + {z} {x [y\ } 

+ {x}{y}{z} (1.2) 

can be extended to higher-order products xo^i . . .x n -i. 

These identities make it possible to answer questions about the distribution mod 1 of 
sequences having the form 

a x n\a 2 n . . . \a k - X n \_OL k n\ \ . . . \ , n = l,2, .... (1.3) 

Such sequences are known to be uniformly distributed mod 1 if the real numbers l,a±, 
. . . , ak are rationally independent [1]; we will prove that (1.3) is uniformly distributed in 
the special case <x\ = «2 = • • • = ®k = ot if and only if a k is irrational, when k is prime. 
(It is interesting to compare this result to analogous properties of the sequence 

a [ain\ [a 2 n\ . . . [akn\ , n = l,2..., (1.4) 

where «o, «i, • • • , «fc are positive real numbers. If k > 3, such sequences are uniformly 
distributed mod 1 if and only if ao is irrational [2].) 
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2. Formulas for the product x xx . . . x n _x. The general expression we will derive for 
xqXi . . .x n -i contains 2 n+1 — n — 2 terms. Given a sequence X = (xo,xi, . . . ,x n -i) we 
regard x n+ j as equivalent to Xj, and for integers a < b we define 

, f 1 , if a = b: 

X a:b =\ ,, 1U| (2.1) 

\x a [X^ +i y- b \ , otherwise. V ' 

Thus X 1:4 = xi [x 2 [x 3 \ J and X 4: ( n+1 ) = x 4 |x 5 [ . . . |_x n _i [^oj J • • • J J • Using this nota- 
tion, we obtain an expression for xoXi . . .x n -\ by taking the sum of 

{X Si:S2 } {X S2:S3 } ... {x Sk: ( Sl+n ^} — ( — l) h [X Si:S2 \ [X S2:ss \ ... [x Sk: ( Sl+n ^\ (2.2) 

over all nonempty subsets £ = {s\, . . . , Sk} of {0, 1, . . . , n — 1}, where si < • • • < Sfc. This 
rule defines 2 n+1 — 2 terms, but in the special case k = 1 the two terms of (2.2) reduce to 

|j^si:(si+n)| _|_ p^si:(si+n)j _ j^si:(si+n) ^2_3^ 

so we can combine them and make the overall formula n terms shorter. The right-hand 
side of (1.2) illustrates this construction when n = 3. 

To prove that the sum of all terms (2.2) equals XoXi . . .x n _i, we replace {X a:b } by 
X a:b — [X a:b \ and expand all products. One of the terms in this expansion is xqXi . . . x n -i; 
it arises only from the set S = {0,1, . . . ,n — 1}. The other terms all contain at least one 
occurrence of the floor operator, and they can be written 

x Ul . . . x Vl —i \_X 1 2 J x U2 . . . x V2 —i \_X 2 3 J x U3 . . . x V3 — i ■ ■ ■ [X ^ J (2-4) 

where u\ < v\ < u 2 < v 2 < 113 < ■ ■ ■ < Vk < n. We want to show that all such terms 
cancel out. For example, some of the terms in the expansion when n = 9 have the form 

Xl [X 2:4 \ x 4 x 5 [X 6:7 \ [X 7 -- 10 \ = Xl [x 2 [x 3 \ J x 4 x 5 [x 6 \ [x 7 [x 8 [x \ J J , 

which is (2.4) with u\ = 1, v\ = 2, u 2 = 4, v 2 = 6, 113 = V3 = 7. It is easy to see that 
this term arises from the expansion of (2.2) only when S is one of the sets {1, 2, 4, 5, 6, 7}, 
{1, 4, 5, 6, 7}, {1, 2, 4, 5, 7}, {1, 4, 5, 7}; in those cases it occurs with the respective signs — , 
+, +, — , so it does indeed cancel out. 

In general, the only sets S leading to the term (2.4) have S = { s \ Uj < s < Vj }U{ Vj | 
uj = v j} UT, where T is a subset of U = { v j | Uj 7^ vj }. IfU is empty, all parts of the term 
(2.4) appear inside floor brackets and this term is cancelled by the second term of (2.2). 
If U contains m > elements, the 2 m choices for S produce 2 m_1 terms with a coefficient 
of +1 and 2 m_1 with a coefficient of —1. This completes the proof. 
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Notice that we used no special properties of the floor function in this argument. The 
same identity holds when [x\ is an arbitrary function, if we define {x} — x — [x\. 
The formulas become simpler, of course, when all equal. Let 

X ~\x[x^\, if > 0; { ] 

and let 

a k = {x-- k }, b k = [x*\. (2.6) 
Then an identity for x n can be read off from the coefficients of z n in the formula 

xz a\z + 2a 2 z 2 + 3a 3 z 3 H , b x z + 2b 2 z 2 + 3b 3 z 3 j 



1 — xz 1 — aiz — a2-2 2 — CI3Z 3 — ■ ■ ■ 1 + 61Z + biz 2 + + 
which can be derived from (2.2) or proved independently as shown below. For example 

x 2 = a\ + 2a 2 - b\ + 2b 2 ; 

x 3 = a\ + 3ciici2 + 3a3 + b\ — 3bib 2 + 363 ; 

x = a 1 + Aa x a 2 + Aa\a 3 + 2a 2 + 4a4 



b\ + Ab{b 2 - 46i6 3 - 2Y 2 + 4b 4 . 



In general we have 



x n =p n (a 1 ,a 2 ,...,a n ) - p n {-b u -b 2 , ■ ■ ■ , -&„) , (2.8) 
where the polynomial 

(ki + k 2 + ■ ■ ■ + k n — 1)1 n h. t, , . 

p n ( ai ,a 2 ,...,a n )= Yl — rrn a k M 2 ---< n (2-9) 

, , „, f— , , Ki! k 2 \ ... k n \ 

ki +2k2~\ \-nk n =n 

contains one term for each partition of n. 

It is interesting to note that (2.7) can be written 

zd 1 zd 1 zd 1 
in = — In — in 



dz 1 — xz dz 1 — a\z — a 2 z 2 — ■ ■ ■ dz 1 + b\z + b 2 z 2 + 
hence we obtain the equivalent identity 

1 _ 1 + hz + b 2 z 2 + b 3 z 3 + ■■■ 
1 — xz 1 — a\z — a 2 z 2 — a 3 z 3 — ■ ■ ■ 



(2.10) 
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This identity is easily proved directly, because it says that + bk = xbk-i for k > 1. 
Therefore it provides an alternative proof of (2.7). It also yields formulas for x n with 
mixed a's and 6's, and with no negative coefficients. For example, 

x 2 = a\ + a 2 + a-ibi + b 2 ; 

x 3 = a\ + 2a\a 2 + a 3 + (a? + a 2 )h + aib 2 + 63 ; 
x 4 = a\ + 3a\a 2 + 2ai<23 + a 2 + a<± + (af + 2a\a 2 + a 3 )&i 
+ (af + a 2 )b 2 + ai& 3 + 64 . 

3. Application to uniform distribution. We can now apply the identities to a 
problem in number theory, as stated in the introduction. Let [0 . . 1) = { x | < x < 1 }. 

Lemma 1. For all positive integers k and I, there is a function fk,i(yi, Vi-, ■ ■ ■ , Vk-i) from 
[0 . . I)*" 1 to [0 . . 1) such that 



x -.k x k 



I kl ikJ \ I A-!/ 



x 



X 

kU 



X 



k-1 



kl I 



(mod 1) 



Proof. Let 



p n (ai, a 2 , • • • , a„_i) = p n {a u a 2 ,...,a n )-n a n 
be the polynomial of (2.9) without its (unique) linear term. Then 
x : ^ x^ 1 1 

~T = ~ki~ ~ki " ■ ' afe_1 ) + ^ fc (~ &1 ' • • • ' -^-1) • 



(3.1) 



(3.2) 



(3.3) 



We proceed by induction on k, defining the constant /i 5 / = for all /. Then if yj = {x J /kl 1} 
and lj = kl l/jl we have 



a. 



X/' 



{lj ((j - ly.yj - f jtl . (yi, . . . , % _i)) } 



and 



bj = 



X' J 

h J 



= u 



xr 



Ij-i 
i=i L 



h J h J 



1,-1 

E 



{ (J - !) ! - /jMj (yi, • • • , } + r 



i=l 

because of the well-known identities 



j J 



(mod kl) , 



«-i 



{lx}={l{x}}, [_lx\=J2l x + i / l \ 



i=0 



(3.4) 



when / is a positive integer. Therefore (3.1) holds with 



fk,i{vu Vk-i) = <| jn Pk{ai,k,i, • • • , a>k-i,k,i) ~ Yi Pk(~ b i,k,h • • • , ~h-i,k,i) j , (3.5) 



where 



aj, k ,i = { (U ~ 1)1 yj ~ fj,kii/ji(yi, • • • , Vj-i))M l/jl} , 

h,k,l 



k\i/j\-i 

E 



i=l 



{(j -!)!%• ~fj,k\i/j\ (?/i> 



(3.6) 
(3.7) 



For example, 



/ 2)3 (?/) = {(«?-/3 1 2 )/6}, 
/ 3 ,i(y, z) = {(30^2 + a? - + ^?)/3} , 

where ot\ = {6z/}, a 2 = {3z - 3 f 2,3(1/)}, Pi = [y + gj + [y + §J 

A = Li> - /2, 3 (y)} + §J + IU - /2,s(y)} + fj. □ 



+ + IJ , and 



Lemma 2. The function fj-j of Lemma 1 does not preserve Lebesgue measure, and neither 
does {klmfk,i} for any positive integer m. 

Proof. It suffices to prove the second statement, for if f^ i were measure-preserving the 
functions {mfk,i} would preserve Lebesgue measure for all positive integers m. Notice that 
{klmfk,i} = {mpk(ai,k,i, ■ ■ ■ ,a,k-i,k,l)}, because Pk(-h,k,l, ■ ■ ■ , -h-i,k,i) is an integer. 
The triangular construction of (3.6) makes it clear that ai,fc,z, • • • , afc-i,fc,z are independent 
random variables defined on the probability space [0 . . l) fc_1 , each uniformly distributed in 
[0 . . 1). Therefore it suffices to prove that {mpk(ai, . . . , at-i)} is not uniformly distributed 
when ai, . . . , cik-i are independent uniform deviates. 
We can express Pk(a>i, ■ ■ ■ , a>k-i) m the form 

k ciia k -i + aiQi(ai, . . . , a fe _ 2 ) + k a 2 a k - 2 + a 2 q 2 (a 2 , • • • , afc-3) H H \k a\ /2 , 

for some polynomials qi, ... , <?L(fc-i)/2j) where the final term |^ a fc/2 * s absent when is 
odd. Then we can let yj = cij for j < ^k and yj = cij — qk-j(cik-j, ■ ■ ■ , a,j-\)/k for j > 
obtaining independent uniform deviates y±, . . . , j/fc_2 for which mpi-( a i, • • • > a fc-i) equals 



</fc(yi, • • -,2/fc-i) = mkyxyk-i +mky 2 yk-2 H h (f ky 2 k/2 [k even]) 



(3. 



For example, ^4(2/1, 2/2, 2/3) = 4yi2/3 + 2j/l and 5-5(2/1, 2/2, J/3, 2/4) = 5yiy 4 + 52/22/3 when m = 1. 

The individual terms of (3.8) are independent, and they have monotone decreasing 
density functions mod 1. (The density function for the probability that {kxy} G [t . . t+dt] 
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is X^j=o \ m <^-) Therefore they cannot possibly yield a uniform distribution. For 
if f(x) is the density function for a random variable on [0 . . 1), we have E(e 27TlX ) = 
J Q e 27TZX f(x) dx y£ when f(x) is monotone; for example, if f(x) is decreasing, the imagi- 

1 /2 

nary part is f Q sin(27rx) (f(x) — /(l — x)) cte > 0. If Y is an independent random variable 
with monotone density, we have E(e 2 ^ x+Y ^) = E(e 27Ti( - x+Y ^) = E(e 27TiX )E(e 277iY ) ^ 0. 
But E(e 27TlU ) = when U is a uniform deviate. Therefore (3.8) cannot be uniform mod 1. 
□ 

Now we can deduce properties of sequences like 

(an) :k = an [an [ . . . [an\ ...J 
as n runs through integer values. 

Theorem. If the powers a 2 , . . . , a fc_1 are irrational, the sequence {m(an) h — km(an) :k }, 
for n = 1, 2, . . . , is not uniformly distributed in [0 . . 1) for any integer m. 

Proof. This result is trivial when k = 1 and obvious when k = 2, since {(an) 2 — 2(an) :2 } = 
{an} 2 . But for large values of k it seems to require a careful analysis. By Lemma 1 we 
have 

{ m (an) k -km(any k } = ^kmf k)1 {^) j - - })} , (3.9) 

and Lemma 2 tells that {kmfk,i} is not measure preserving. 

Let S be an interval of [0 . . 1), and T its inverse image in [0 . . l)^ -1 under {kfk,i}, 
where /z(T) 7^ A*(5')- It is easy to see that if (yi, . . . , j/fc-i) £ ^ and yi,...,yjfc-i are 
irrational, there are values ei, . . . , e^-i such that [yi . . yi + ei) x • • • x [y/c-i • • Vk-i + 
efc-i) ^ T. Therefore the irrational points of T can be covered by disjoint half-open 
hyperrectangles. We will show that (3.9) is not uniform by using Theorem 6.4 of [3], 
which implies that the sequence {{a\n ei },..., {a s n es }) is uniformly distributed in [0 . . 
l) s whenever a±, . . . , a 8 are irrational numbers and the integer exponents ei, . . . , e 8 are 
distinct. Thus the probability that {(an) k — k(an) :k } G S approaches /z(T) as n — > 00; 
the distribution is nonuniform. [] 

Corollary. If the powers a 2 ,...,a fc_1 are irrational, the sequence {{an) :k }, for n = 
1,2,..., is uniformly distributed in [0 . . 1) if and only if a k is irrational. 

Proof. If a k is irrational, {a k n k /k} is uniformly distributed in [0 . . 1) and independent 
of ({an/k\}, . . . , {a k ~ 1 n k ~ 1 /k\}), by the theorem quoted above from [3]. Therefore the 
right-hand side of (3.1) is uniform. 

If a k is rational, say a k = p/q, assume that {{an) :k } is uniform. Then {q(a k n k — 
k(an) :k )} = {—qk(an) :k } is also uniform, contradicting what we proved. □ 
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We conjecture that the theorem and its corollary remain true for all real a, without 
the hypothesis that a 2 , . . . , are irrational. 
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